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Components for Optical Qubits in the Radio Frequency Basis 
5 Technical Field 

The invention is in the field of optical quantum information processing. This invention 
concerns the manipulation of single photon qubits in frequency modes, and in particular 
a frequency beamsplitter and a frequency half-wave plate for use in the frequency basis. 

10 

Background Art 

Quantum information operations may be performed by assigning logical states to the 
quantum states in physical systems. The logical states are analogous to the zero and 

15 one states for a conventional digital computer. Quantum information can be encoded 
and manipulated using single photon states. Many in principle demonstrations of 
quantum information tasks have now been accomplished in single photon optics 
including quantum cryptography [1], quantum dense coding [2] and quantum 
telepprtation [3]. More recently two qubit gates have been realised [4,5] using 

20 conditional techniques [6]. Such experiments typically make use of polarisation to 
encode the qubits. 

The tools used in polarisation encoding schemes are principally the half-wave plate 
which is used to make arbitrary rotations of a state in the polarisation basis, the 
25 polarising beamsplitter which is used to separate (or combine) photons into (from) 
different spatial and polarisation modes and the quarter-wave plate which is used to 
introduce relative phase shifts between the two bases. 

It will be appreciated that polarisation is not the only photonic degree of freedom 
30 available to the experimentalist. For example schemes in which the temporal [7] or 
spatial [8] optical modes form the quantum basis have also been realised. Here we 
consider a basis consisting of two frequency optical modes. 



35 



Pi - » 

WO 2005/043787 



PCT/AU2004/001476 



2 . 

Disclosure of the Invention 

In a first aspect the invention is a frequency beamsplitter for use in the frequency basis,, 
comprising: An asymmetric two-path interferometer, reversible down to the quantum 
5 limit. A first partially transmitting mirror to split photons into first and second paths. A 
time delay element to introduce a differential time delay into the second path. And a 
second partially transmitting mirror to mix the two paths again to form two outputs. 

Such a frequency beamsplitter separates two frequency modes on a single optical beam 
1 0 into two spatially separated beams using standard optical elements. 

Suitable interferometers include the Mach-Zehnder interferometer, a Michelson 
interferometer or any other interferometer with a frequency dependent phase shift 

1 5 An optical qubit can be encoded in the occupation by a single photon of one of two 
different frequency modes These frequencies would ideally be sufficiently close 
together that they can be manipulated with standard electro-optical devices but still be 
clearly resolvable using narrowband optical and opto-electronic systems. 

20 The beam splitter may also be implemented using commercial optical fibre or planar 
waveguide technologies. Hence, such an encoding scheme is attractive from the 
perspective of developing stable, robust and ultimately commercially viable optical 
quantum information systems. 

25 The attraction of optical fibre based quantum optical systems is clear. For example, 
there is ongoing interest in developing non-classical optical sources that will be well 
suited to optical fibres and optical fibre technologies [9-11]. Indeed a Quantum Key 
Distribution (QKD) scheme using radio frequency amplitude and phase modulation as 
the conjugate bases has recently been demonstrated using optical fibres and fibre 

30 technologies [12]. 

In free space, additional mirrors may be provided between the partially transmitting 
mirrors or at the output to redirect the first and second paths to provide the outputs at 
convenient locations. Optical fibre or planar waveguide structures may perform similar 
35 functions. 
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In a second aspect the invention is a radio frequency half-wave plate for use in the 
frequency basis, comprising: A frequency beamsplitter as defined above to split a 
beam of photonic qubits encoded in different frequency modes. An acousto-optic 
modulator driven to receive the outputs from the frequency beamsplitter and couple the 
5 two different frequency and spatial modes of the outputs together via a phonon 
interaction to produce two further outputs. A frequency beamsplitter as defined above 
operating in reverse to receive the two further outputs from the acousto-optic. modulator 
and produce a combined output. 

10 Such a device is able to produce any selected rotation of photonic qubits in the radio 
frequency basis. 

It should be appreciated that in the frequency basis relative phase shifts, that is quarter 
wave plate (QWP) action, is naturally achieved during propagation. 

15 

In summary, the invention provides devices which may be used as the principle 
experimental components in optical quantum information systems which make use of 
the frequency basis. These components are essentially the frequency basis analogues of 
the polarising beamsplitter and the half-wave plate. They allow arbitrary, reversible 
20 operations to be perfomed on frequency basis encoded qubits. 

Brief Description of the Drawings 

An example of the invention will now be described with reference to the accompanying 
25 drawings, in which: 

Fig. 1(a) is a graph illustrating | l) z , Fig. 1(b) is a graph illustrating 1 0) z and 
Fig. 1(c) is a graph of the sum of 1 l)^ and | 0) z , namely 

Fig. 2 is a schematic diagram of a radio frequency beamsplitter. 
Fig. 3 is a . schematic diagram of a radio frequency half-wave plate incorporating 
30 the beamsplitter of Fig. 1 . 

Fig. 4 is a folded form of the radio frequency half-wave plate of Fig. 3. 
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Best Modes of the Invention 

Referring first to Fig. 1(a) and (b) to illustrate the physics of the invention, we first 
present an idealized example of an encoding scheme which makes use of the radio 
5 frequency basis. Consider the logical basis whereby 



10 where | 0) L and [1)^ denote the logical states of the qubit and the. notation \n) 

denotes an n photon state at the frequency ±Q. relative to the average or carrier 
frequency a> 0 . We shall assume that the states are indistinguishable apart from their 

frequencies. Given current technical considerations, CI is taken to be a radio or 
microwave frequency in the range of tens of megahertz to a few gigahertz. It is 
1 5 convenient to write the states of Eq.l in the form 

|0) X =^(-Q)*]0) 

|l),=^4G)t|0) (2) 

where A^co)* is the creation operator for the a> frequency mode. As all the elements in 
20 our device are passive (energy conserving) we can obtain the state evolution produced 
through the device by considering the Heisenberg evolution of the relevant operators. 

Fig. 2 is a schematic diagram of a radio frequency beamsplitter comprising a highly 
asymmetric Mach-Zehnder interferometer 20. This is arranged with an input port 22 to 

25 receive a radio frequency input signal A^. The arrows indicate directions of 
propagation of internal optical fields. A first partially transmitting mirror 24 which 
splits the optical input A^ into first 26 and second 28 paths. First 30 and second 32 
mirrors redirect the first 26 and second 28 paths. In the second path 28 there is a time 
delay element 34 to introduce a differential time delay r into the second path. The two 

30 paths are then mixed again at a second partially transmitting mirror 36 to form two 
outputs A x and A 2 . 
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The annihilation operators for the forward traveling outputs of the Mach-Zehnder 
interferometer at the frequency co may be written as 

4 (a) = |[4, (pX-e^e™ + 1) + {a>){e^e w + 1)] 

! (3) 
^ (fl>) = -[^ (*>)(>V~ + 1) + » fc - 1)] 

where ^ and ^ are defined in Fig. 1, r is the differential time delay introduced into 
one of the arms of the interferometer and <j> = a> 0 r is the phase shift acquired by a field 
at the carrier frequency. 

10 Choosing the time delay r in the interferometer such that (j) =7r/2 and Cir = *z72, the 
creation operators for the forward travelling outputs of the Mach-Zehnder 
interferometer at the frequencies co = ±Q are 

ACQ)* =^(£2)', A 1 (^f2) t =-io. 7l (-n) t 

15 

In the state picture we have that an arbitrary input state 

|#>-(M,(-0)t +v^(Q) t )|0)^|0) u;n (4) 
20 is transformed to the output state | ip') according to 



>') = u\ ip) = t/K (-«)* + vA, («) + )l 0) ^ 1 0) 

= O/tt^C-n)^* +vc/4„(^)V t )l o) j o) 
= (MiA 2 (-ay +vA l (ay)\o) Al \o} 



A2 

A2 

A2 



(5) 



where U is the unitary operator representing the evolution through the element. In 
25 going from lines two to three we have used the fact that UA bt (-O) 1 " Z7 f is time reversed 

Heisenberg evolution, obtained explicitly by inverting the standard Heisenberg 
equations such that the input operator is written in terms of the output operators. We 



WO 2005/043787 



6 



PCT/AU2004/001476 



have also used t/j 0)^10)^ =\0) Al \0) A2 . The action on the frequency encoding is 
equivalent to the action of a polarising beamsplitter in polarisation encoding. 

The frequency beamsplitter 20 is used as the first stage of the radio frequency half- 
5 wave plate (RF-HWP) 40 which will now be described with reference to Fig. 3 . 

The radio frequency half^wave plate (RF-HWP) 40 of Fig. 3 comprises a first frequency 
beamsplitter 20 followed by an acousto-optic modulator (AOM) 42, followed by 
another frequency beamsplitter 20\ The arrow in the AOM 42 indicates the direction of 
10 propagation of the acoustic wave in the device. The frequency of the radio frequency 
source 44 used to drive the AOM is S . 

In the Heisenberg picture, we define the annihilation operator for the input mode at a 
particular Fourier frequency co relative to. the carrier as A^ipo) . We also define an 
1 5 ancilla field u fa (t) , initially in a vacuum state, entering the device vertically from the . 

bottom. Fig. 3 defines a number of internal annihilation operators for the RF-HWP as 
well as two output fields. The output of interest is A m . 

It can be seen that, in use, output A, is reflected from a further mirror 46 to a first input 
20 port of the AOM, and output A 2 also passes into to the second input port of the AOM. 

If we take the logical basis as defined in Eq. 1 and 2, we become interested specifically 
in the operators, A k (Q) and A k (-Cl) where k = in,l..6,out,back . We shall make use 
of the relation that Al(w) = A k (-co) f [16] to find the relevant creation operators. 

25 

The heart of the RF-HWP is the acousto-optic modulator (AOM) 42 which couples two 
different frequency and spatial modes together via a phonon interaction [17,18]. In our 
scheme the AOM is used to shift photons between the two frequencies -Q and f2. 
The an n i hil ation operators for the outputs of the AOM, A 3 (co) and A A (co) as defined in 

30 Fig. 3, are [17] 

A 3 (co ) = cos 0A X (co ) + 1 sin 0A 2 (co - S) (6) 
A A (co) = cos 0A 2 (co) + z sin 6A X (co + S) (7) 

35 where S represents the modulation frequency applied to the AOM and 0 is a measure 
of the diffraction efficiency of the AOM such that cos# represents the undiffracted 
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fraction of the field and sin£ represents the diffracted fraction.. We have taken the 
asymmetric phase convention for the AOM outputs and note that 6 is proportional to 
the amplitude of the radio frequency modulation applied to the AOM [18]. We note that 
the second term in Eq. 6 represents the diffracted, and hence frequency upshifted, 
5 component of the input field A^. Similarly, the second term in Eq. 7 represents the 
downshifted component of A 1 . 

The output fields A 3 , after reflection in mirror 48, and A A pass through the second 

Mach-Zehnder interferometer 20'. The annihilation operators for the output fields of the 
10 system, A oul and A back , are 

A M (<o) = ^[A 4 (aX-e* 2 e™*+l) + iA 3 ((Q)(e i **e™*+l) (8) 
A l ^(p) = ^\iAi(a>)(e**e vm * +l)+A 3 (oi>Xe l *V ,ot * -1) (9) 

1 5 where <j> 2 = an 2 is the phase shift acquired by the "carrier" in the second interferometer 
with time delay x 2 . We can combine the results of Equations 3, 6 and 7, 8 and 9 as 
well as making use of A* (o) = A(~Q)) 1! to arrive at the creation operators for the 
outputs of the RF-HWP in terms of the inputs. Focussing on the output A^ at the 

frequencies of interest co = ±Ci, setting the AOM modulation frequency S = 2£l and 
20 setting the time delays in both interferometers such that ^=tt/2, (J> 2 =7c/2 and 
Qx — Q.x 2 —n/2, we find that 

A^C-Ciy =-r[cos9A in (-n) t + sin0A in (Q) t 
A out = -i[cos9A,„ (Q) f - sin0A M (-CI? 

25 Applying Eq. 12 to the | 0) L input state leads to 

^>o. ow =-[cosBA oul (-Ciy -sineA^^JO^IO) 
= -(cosei0) i -sineil> i .) 



30 



(10) 



DOUt 



(ID 



whilst applying it to the 1 l) £ input gives 

= -(COSO 1 1)^ + 8^1910),.) 
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Equations 11 and .12 are the key results. Up to a global phase, these equations are 
formally equivalent to those used to describe the rotation of an arbitrary two 
dimensional vector through the angle 0. Hence the system illustrated in Fig. 3 is 
5 operationally equivalent to a half-wave plate in the basis defined by the frequencies 
- Q, and Q. . 

Let us now turn to technical issues. Conceptually, the RF-HWP comprises a FBS, 
followed by an AOM followed by another FBS. However, implementing the RF-HWP 
10 in such a fashion would be a tremendously challenging technical task. It would require 
actively locking the phase 0=n/2 in two different interferometers. In addition, the 

optical path length between the two interferometers would need to be locked. This is 
why we have proposed the RF-HWP with the folded design illustrated in Fig. 4. This 
folded design requires locking of only one interferometer. Further, a locking signal can 
15 be derived from, the backward travelling output of the RF-HWP, A badc , without 
disturbing the useful output A aa 

Let us denote the "forward travelling" beams as those propagating to the right of Fig. 1 
and the "backward travelling" beams as those propagating to the left. 

20 

We double pass the AOM in this scheme by reflecting A 3 in mirror 50 and A4 in mirror 
52 . The two backward travelling fields A s and A 6 emerge from the AOM as illustrated 

in Fig. 5. The annihilation operators for the backward travelling fields emerging from 
the AOM are 

25 

A 5 (a) = cos0A A (to) + 1 sin0A 3 (a + S) 

= cos 20A 2 (w) + 1 sin 20A X (co + S) ^ 
A 6 (a)) = cos0A 3 (o) + isia0A 4 ((v-S) 

= cos 20A X (w) + 1 sin 20A 2 (co - S) ^ 

The backward travelling fields A 5 and A 6 make a second pass of the Mach-Zehnder 
30 interferometer 20. The annihilation operators for the output fields of the system, A oja 
A back> are 

A owf (cD) = ^[iA 5 (cD)(^^ +l) + A 6 (a>)(eV- -1) (15) 
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^(^) = |k(«)(-^e'' Br +l) + ^ 6 (a>)( e '< i e ^+l)] (16) 

where the parameters <j> and t are the same as those for the forward travelling 
asymmetric Mach-Zehnder interferometer. 

5 

We can combine the results of Equations 3, 13, 14, 15 and 16 as well as m aking use of 
A*(g>) = A(-(oy to arrive at the creation operators for the outputs of the RF-HWP in 
terms of the inputs. Focussing on the downward travelling output A m at the 

frequencies of interest co = +Q., setting the AOM modulation frequency 5 = 2Q and 
10 setting t such that <j>-nl1 and Clr = n 1 2, we find that 



15 



20 



A oia (-ay =-[cos29A &l (-a) t +sm2QA in (a? 
A^iCl? =-[cos28A i> ,(n) t -sin2eA £ „(-a) t 

Applying Eq. 19 to the ) 0) z input state leads to 

I *) 0 ,o,a = -[cos2^(-n)T -sin^CQ)^! 0)^1 0} 
= -(cos2<9| 0) L -sin l) L ) 

whilst applying it to the 1 1) £ input gives 

I *)o, out = -[cos2M 0 « (^) + +sin20A mu (-ay]\ 0) A J 0) 
= -(cos 20\ 1) L + sin 20\ 0) L ) 



(17) 



uout 



(18) 



VOtlt 



(19) 



Equations 18 and 19 are again the key results. Up to a global phase, these equations are 
formally equivalent to those used to describe the rotation of an arbitrary two 
dimensional vector through the angle 0 = 20 . Hence the system illustrated in Fig. 5 is 
25 operationally equivalent to a half-wave plate in the basis defined by the frequencies 
-CI and Q. 



30 



The situation considered so far is of course impractical, as single frequency qubits will 
be stationary in time. More realistically, we might consider finite bandwidth qubits of 
the form: 
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I <>>x « (— > V4 l^e-^^A in (coy\ 0) 

™ (22) 

The overlap between these qubits is 
5 (0/l) L =6- 2 ° 2 ^ (23) 

thus provided the width of the frequency packet is sufficiently small compared to its 
mean (i.e a « 2Q. 2 ) then these qubits will be approximately orthogonal. The problem 
with the finite frequency spread for our device is that now the condition Qr = n 1 2 will 
10 not be precisely satisfied for the entire frequency packet. The effect is to produce a 
phase shift across the frequency wave packet and also to produce some probability of 
photons exiting the device in the wrong beam (ie A back ) or at frequencies outside the 

computational basis (eg 2Q). Taking such events as lost photons and tracing ovqr them 
leads to a mixed state which can be written 

15 



20 



25 



Po>*j=\Q)AQ\i+\Q)\Q\ (?4) 

where , represents the mixed state output obtained for the logical state input | i) L . 
For the 1 0) x input state 

| <2) n = -(— ) 1/4 f ^e- (n+a,)I/flr (cos2^4„(G?) t l/2(l + e" r<n+a,)/n )- 

tct 1 ^ (25) 

smWA^ (co + 2Q) + 1/ 4(1 + e '" n ^' n f )| o) 
and for the 1 1) £ input 



_2_ 
na 

sin 2fl4 ftl (o - 2<Q) t l/4(l + e -**ivp~v*^ 0 ) 



| Q\ = -(^r) 1M J da>e-^- af ' a (cos 20^ (o>) t 1 / 2(1 + e ->*< n -«)'° ) + 



and where \Q} is a collective ket representing all the photons that end up in 
(orthogonal) states outside the computational basis. We can evaluate the impact of this 



WO 2005/043787 



PCT/AU2004/001476 



11 

effect by calculating the fidelity between the expected output state, \P) }i and that 
obtained: 

F-(*\'P~j\*) § (27) 

5 

The expression for the fidelity is complicated and of limited utility to reproduce here, 
however its pertinent features can be listed succinctly: the fidelity depends strongly on 
the ratio of CI 2 to or ; it depends weakly on the rotation angle 0 ; as expected it tends to 
one as the ratio Q. 2 la tends to infinity. Some representative results are: Q 2 la = 10 , 
10 F = 0.942; Q 2 /<x = 100, F = 0.9939; Q 2 /a =1000, F = 0.99938 . We conclude that 
high fidelities are consistent with sensible signal bandwidths. 

The technical limitations to the performance of the RF-HWP will be set by the 
diffraction efficiency of the. AOM, the transmission losses of the AOM and the mode- 
ls matching efficiency in the interferometer. Rotation of the input through an angle of 
© = nil requires that the AOM have a diffraction efficiency of 50%. This technical 
requirement can currently be met with commercially available devices [21]. 
Transmission losses in the AOM mean that the mode which is redirected into the 
interferometer, and which is ultimately measured, is actually a mixture of the original 
20 mode and the vacuum. This effect may be treated as a perfectly transmitting AOM 
with a partially transmitting beamsplitter placed on each output port [19,20] Mode 
mismatch in the interferometer will reduce the fringe visibility at the sideband 
frequencies. . We can model this effect by decomposing all fields into projections onto 
an orthonormal basis set of spatial modes, and then introducing a rotation of one of the 
25 internal fields of the interferometer relative to the other in this vector space. We treat 
photons which appear in the incorrect spatial modes as lost. 

We can evaluate the impact of AOM losses and mode-matching efficiency by 
calculating the fidelity between the expected output state and that actually obtained at 
30 the output of the RF-HWP. The fidelity after taking account of AOM losses and 
mode-matching efficiency, F'is 

p. m i?(W[2(l + Vmm )cos 2 (20)+ (l + sin 2 (2*)] J, (28) 
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where the single-pass power transmission in the AOM is given by tj acm and the 
single-pass mode-matching efficiency of the interferometer is given by tj^. The 

fidelity of the system in the absence of technical limitations is given by F. A good 
quality free-space AOM would have, rj > 0.95 [21]. Similarly, a well mode-matched 
5 interferometer would have rj nm > 035 [22]. We note in passing that detector 

efficiency or homodyne mode-matching efficiency may be modeled in the same way as 
losses in the AOM . 

Suitable interferometers include the Mach-Zehnder interferometer, a Michelson 
1 0 interferometer or any other interferometer with a frequency dependant phase shift. 

It will be appreciated by persons skilled in the art that numerous variations 
and/or modifications may be made to the invention as shown in the specific 
embodiments without departing from the spirit or scope of the invention as broadly 
15 described. The present embodiments are, therefore, to be considered in all respects as 
illustrative and not restrictive. 
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